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We report a detailed investigation on the lower critical held H c \ of the superconducting 
Bao.6Ko.4Fe2As2 (FeAs-122) single crystals. A pronounced kink is observed on the H C \{T) curve, 
which is attributed to the existence of two superconducting gaps. By fitting the data H c i (T) to the 
two-gap BCS model in full temperature region, a small gap of A a (0) = 2.0 ± 0.3 meV and a large 
gap of Af,(0) = 8.9 ± 0.4 meV are obtained. The in-plane penetration depth A a t(0) is estimated to 
be 105 nm corresponding to a rather large superfluid density, which points to the breakdown of the 
Uemura plot in FeAs-122 superconductors. 
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One of the crucial issues in understanding the super- 
conducting mechanism in recently discovered FeAs-based 
layered superconductors [J, Q is the pairing symmetry of 
the superconducting gap and the nature of the low en- 
ergy excitations. For LnFeAsOi-zF^ (FeAs-1 111) super- 
conductors with x = 0.04 ~ 0.2 and Ln=La,Ce,Nd,Sm, 
reports on pairing symmetry are divided into two cat- 
egories: those favoring a gap with 0, 0, 0, 0, 0] or 
without [1, 0, [n| nodes. In other aspect, the upper 
critical field measurement (llj . tunneling spectroscopy 
O, NMR_Jl3|], magnetic penetration depth measure- 



ment 



l4 1 1 51 ] and Mossbauer [16| experiments indicate 
the existence of a two-gap nature. However, because 
most of the current experiments were performed on poly- 
crystalline samples, the experimental results within the 
context of pairing symmetry have not yet reached a con- 
sensus. 

The emergence of (Ba,Sr)i_ x K x Pe 2 As 2 (FeAs-122) 
171 Il8j hole-doped superconductors with transition tem- 
perature T c up to 38 K has enriched the research in 
this area. Although stoichiometric (undoped) BaFe2As2 
shares the same feature with undoped LnFeAsO of a spin 
density wave type magnetic order and a structural transi- 
tion at similar temperatures , however, there are 
also some significant differences between the properties 
of FeAs-1111 and FeAs-122 superconductors. Most no- 
table is the fact that FeAs-1111 superconductors seem to 
have very low charge carriers density Q and hence low 
superfluid density [a, l2l( . This may give explanation to 



the fact that the Uemura plot is satisfied in the FeAs- 
1111 systems [Il|,[22j]. According to Uemura et.al. [23l |. 
in superconductors with low superfluid density (p s ), T c 
scales linearly with p s cx A~ b 2 (0), this is actually not re- 
quired by the BCS theory. However, it is found that the 
charge carrier density in FeAs-122 is an order of mag- 
nitude larger than those in FeAs-1111 systems 171 124|. 
Thus it is very interesting to know whether the Uemura 
plot is still satisfied in FeAs-122 superconductors. 

Lower critical field H c i(T), or equivalently, magnetic 
penetration depth A(T) are fundamental probes of the 



nature of the pairing symmetry, or/and multigap of un- 
conventional superconductors. As an advantage, H c i(T) 
measurement probes relatively large depths (in scale of 
A ~ 100 nm) and are insensitive to sample surface 
condition. In this Letter we present the first detailed 
magnetic penetration measurements of superconducting 
Bao.6Ko.4Fe2As2 single crystals. The local magnetiza- 
tion measurements allow a precise determination of H c \ . 
We found the presence of possibly a full gap feature to- 
gether with two gaps in Bao.6Ko.4Fe 2 As2 superconduc- 
tors. Meanwhile, the absolute value of H c i(0) determined 
from this work place the samples far away from the Ue- 
mura plot. 

Crystals of Bao.6Ko.4Fe2As2 were grown by FeAs flux 
method. The details of crystal synthesis have been de- 
scribed elsewhere 24|. Our crystals were characterized 
by both resistivity and ac susceptibility measurements 
with T c — 36.2 K and a transition width of AT C = 0.45 
K (10%-90% of normal state resistivity) [13], as displayed 
in left inset of Fig. 1 (a). A single crystal (sample No. 1) 
was selected from the cleaved as-grown bulk under opti- 
cal microscope. The sample has dimensions of 110 /im 
in diameter and 40 /im in thickness, as shown in the 
right inset of Fig. 1(a). The crystal structure was exam- 
ined by X-ray Diffraction (XRD), and a typical diffrac- 
tion pattern was shown in the main panel of Fig. 1(a). 
Only (00/) peaks were observed with the full-width-at- 
half- maximum (FWHM) around 0.1 °, which indicates 
good crystallization of our samples. 

The local magnetization measurement was performed 
on two crystals using a two dimensional electron gas 
(2DEG)-based micro Hall sensor with an active area of 
10 x 10 /im 2 . The Hall sensor was characterized without 
sample attachment at different temperatures, the Hall 
coefficient Rr was identical to 0.22 il/Oe, independent 
of T below 80 K. In our experiment, we used a small field 
sweep rate of 30 Oe/min to measure the isothermal mag- 
netization M(H) curves in both decreasing (Md ec ) and 
increasing (Mi nc ) the fields to minimize the complex ef- 
fects of the character of the field penetration in a layered 
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FIG. 1: (Color online) (a) X-ray diffraction pattern of the 
Bao.6Ko.4Fe2As2 crystal. The right inset shows the sample 
(No. 1) under the optical microscope. The scale underneath 
with 100 ^im/grid is used to display the size of the sample. 
The left inset shows the resistive transition on one crystal 
cut from the same piece as the one for the measurement of 
H c i- (b) The induction B by Hall probe as a function of T 
measured in zero-field-cooled mode with applied field. Inset: 
A typical magnetization hysteresis loop measured at 32.1 K 
by the Hall sensor. 



structure [25|. A second crystal (sample No. 2) with di- 
mensions of 210 x 150 x 50 /im 3 was also measured by Hall 
sensor and vibrating sample magnetometry (VSM), and 
both measurements showed essentially identical behavior 
except for the different demagnetization effect. 

In order to check whether the vortex entry is influenced 
by the Bean-Levingston surface barrier (2rj 27|. we mea- 
sured the magnetization hysteresis loops. The result is 
shown in the inset of Fig. 1(b). The M(H) curve shows 



a symmetric feature at T = 32.1 K which is close to T c , 
indicative of the dominance of bulk pinning and the ab- 
sence of the Bean-Levingston surface barrier for vortex 
entry. 

The quality of the crystal and the accuracy of the local 
magnetization measurement were demonstrated in Fig. 
1(b), in which the induction B = H + 4ttM sensed by 
the Hall probe is plotted as a function of T under a bias 
field H = 20 Oe. The superconducting transition can be 
detected as the induction jump at T c — 35.8 K with a 
width of AT C = 0.5 K, which shows the high quality of 
the crystal. As shown in raw data of B(T), at low tem- 
peratures, B is close to 0, implying a full Meissner shield- 
ing effect: AirM ~ —H, Alternatively, the magnitude of 
the change in B measured by the Hall sensor through 
the superconducting transition is 19.8 Oe, very close to 
the applied magnetic field 20 Oe. Thus the achievement 
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FIG. 2: (Color online) (a) A typical magnetization hystere- 
sis loop (symbols), the solid lines are the linear fitting curves 
using the low field data (Meissner line). The increasing and 
decreasing branches are shifted downward and upward, re- 
spectively for clarity. The same case is shown in the bottom- 
left inset when the maximum field is less than H c %. The 
upper-right inset shows the magnetization data subtracted 
the Meissner line. The arrows indicate the direction of sweep- 
ing fields and the determination of H&. (b) The initial part 
of the magnetization curves M(H) of sample No. 2 at vari- 
ous temperatures. The dashed line gives the Meissner linear 
approach. Inset: The same magnetization data in field sub- 
tracted by the Meissner line. The dashed line in the inset sets 
up a criterion of 0.5 Oe. 



of full Meissner shielding effect in our measurement pro- 
vides a reliable way to determine the value of H c \ . 

Shown in Fig. 2(a) are the typical isothermal M(H) 
curves by taking Mdec and Mi nc at T=17.1 K, respec- 
tively. It can be seen that, at low H, the Mdec and Mi nc 
arc fully reversible, showing a common linear dependence 
of the magnetization on field as displayed in the bottom- 
left inset of Fig. 2(a). At high H, a deviation from 
the linear dependence occurs at H = H c i for both Md ec 
and M inc curves. To distinguish these deviations, we fit 
more than 50 data points between 10 and 30 Oe by a 
linear relation to account for the common linear depen- 
dence of M(H). These fitted linear lines describe the 
Meissner shielding effects ("Meissner line") at low fields, 
as evidenced quantitatively in Fig. 2(a) in which the 
slope of the fitted lines are -0.98, being very close to - 
1. Thus the deviation of M(H) from the linear Meissner 
shielding is an indication of the first penetration field 
H c i. An alternative way to determine the value of H& 
from these reversible isothermal M(H) curves is to sub- 
tract the Meissner line from both Md ec (H) and Mi nc {H) 
curves, as illustrated in the upper-right inset of Fig. 2(a). 
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The threshold field of non-zero magnetization happens to 
be the divergence of the increasing and decreasing M (H) 
curves. It is noted that the values of H cl with this crite- 
rion were determined in both increasing and decreasing 
field, so they are the true thermodynamic values and are 
not altered by the surface barrier [25 1. 

For a strict treatment, we determined the value of H c \ 
by examining the point of departure from the Meissner 
line on the initial slope of the M(H) curve. In the inset 
of Fig. 2(b) we show how to determine H c i by using a 
criterion of B = AM = 0.5 Oe at different temperatures. 
For a quantum flux O = 20.7 Oe pm 2 , AM = 0.5 Oe 
is equivalent to about (2 ~ 3)00 penetrating into the 
(10 x 10) pm 2 sensing area, which is the limit of our 
Hall probe technique. The H c i values determined in this 
way are about 4% larger than those estimated from the 
point where the reversible magnetization deviates from 
linearity, and we did not observe any significant difference 
in the T-dependence of H c i deduced from either of the 
two criterion. 

Shown in Fig. 3(a) and (b) are the main results of our 
experiment, in which the obtained H c \ are plotted as a 
function of T for crystals No. 1 and 2, respectively. At 
T < 4 ~ 5 K, H c i(T) is weakly T-dependent and seems 
to show a tendency towards saturation at lower T (in the 
limited temperature range). As illustrated in the insets 
of Fig. 3(a) and (b), the saturated H c \ reach 695 Oe for 
sample No. 1 and 590 Oe for sample No. 2. This ten- 
dency of H c i(T) reflects a possible fully gapped nature 
of superconducting state at low T for Bao.6Ko.4Fe2As2 
superconductors, although we could not rule out the pos- 
sibility of a small gap with nodes in the dirty limit. 

A pronounced kink can be easily observed in H C \(T) 
curves at T ~ 15 K for both samples. Obviously, the oc- 
currence of the kink in H c i(T) can not be explained by 
the model with an s-wave or <i-wave single gap. On the 
other hand, this kinky structure in H c \ (T) resembles that 
of the related penetration depth of the two-band super- 
conductor MgB 2 [28|, in which a positive curvature was 
observed and explained by the multi-band theory [29] ]. 
In addition, recent ARPES measurement resolved a two- 
gap nature in a similar Bao.6Ko.4Fe2As2 crystal [30I. l3lj|. 
Thus our observation of a kink in H c i(T) strongly sug- 
gests the existence of multiple gaps in Bao.6Ko.4Fe2As2 
superconductors, being consistent with that predicted in 
electronic band structure calculations [32j ]. 

Quantitatively we analyze and discuss in more detail 
the data of H c i(T) to elucidate the related gap func- 
tion. For our crystals, assuming A a &(0) ~ 100 — 200 nm 
[111], the coherence length £ a b(0) was estimated to be 2- 
2.5 nm from an extremely high upper critical field H C 2 

<A>/2<yo)] 



341 ] . and the mean free 
path (determined from the resistivity at 38 K) is ~ 15 
nm, our samples are therefore expected to be in the mod- 
erately clean, local limit. In this case the local London 
model is valid to describe the data. For a single gap su- 
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FIG. 3: (Color online) The extracted H c \ as a function of T 
for sample No.l (a) and No. 2 (b)(full circles with error bars). 
The open dotted circles in (b) represent the H c \ determined 
from the VSM measurements taken into account of the de- 
magnetization effect. The solid lines are the fitting curves 
using the two-gap model [Eq.(l) and Eq.(2)]. The contribu- 
tions of the small gap [H c i(A a )] and the large gap [H c i(At)] 
in the model are also shown by the dotted and dashed lines, 
respectively. The two insets in (a) and (b) show the enlarged 
views of the data H c i (T) in low temperature region together 
with the theoretical fitting curves (solid lines) with a two-gap 
model. 



perconductor, H c \ relates the normalized supcrfluid den- 
sity as: p,(T) = ^(OVA^r) = H cl (T)/H cl (0), and 



p s (T) is given by 
Ps(T) -- 



35! 3G] 



df(E) 



E 



dE _ A (T) : 



-AE (1) 



with / the Fermi function. Here the total energy is 
E = Ve 2 + A 2 , and e is the single-particle energy mea- 
sured from the Fermi surface. It is assumed that the gap 
A on each Fermi surface follows the weak-coupling BCS 
temperature dependence. For a superconductor with two 
gaps, the normalized superfluid density may be written 
as 



p s = xp s a + (1 - x)p s 



(2) 



where x is the fraction of superfluid density p s a associ- 
ated with the small gap A a . The results of the calcula- 
tions and H c \ of sample No. 1 and 2 are shown by the red 
solid lines in Fig. 3(a) and (b), respectively. Fitting the 
data to above equations (with two s-wave gaps) yields: 
A a = 1.6 ±0.3 meV, A fc = 9.1 ±0.3 meV and x=0.72 for 
sample No.l, and A a = 2.2±0.2 meV, A h = 8.8±0.3 meV 
and x=0.70 for sample No. 2. The gaps obtained from our 
H c i(T) measurements are clearly smaller than those de- 
termined from the ARPES measurements diHH. tws 
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FIG. 4: (Color online) The correlations between T c and the 
superfluid density n a /m* . Points for the cuprates are taken 
from Ref.21,22. The full triangles are showing the data for 
F-LaFeAsO (Ref.5 and Ref.21) and F-SmFeAsO (Ref.13 and 
our work), both are close to the Uemura plot of the cuprate 
superconductors. The data of the FeAs-122 phase (the red 
full circle), MgB2 and NbSe2 are far away from the general 
Uemura plot [231 ]. 



critical field H c i(T) is reliably extracted. It is found 
that H c i exhibits a pronounced kink at T ~ 15 K, 
which indicates a multi-gap nature. By using the two-gap 
weak coupling BCS model to fit the data, we obtained a 
small gap of A a (0) ~ 2.0 ± 0.3 meV and a large gap of 
Af,(0) ~ 8.9 ± 0.4 meV. An estimate of the in-plane pen- 
etration depth gives A a &(0) ~ 105 nm, which points to 
the breakdown of the Uemura relation for the optimally 
doped Bao.6Ko.4Fe2As2 superconductors. 
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discrepancy may be induced by the different ways and 
different criterions in determining the gaps, this should 
be checked by future experiments. It is interesting to note 
that the large gap accounts for only 30% of the total su- 
perfluid density. We must stress that although a small 
gap with nodes (in the dirty limit) cannot be excluded 
from our low temperature data, this will not lead to a 
significant change to the general fitting results obtained 
here. 

Another important result of our experiment, the abso- 
lute value of A a b(0) is independently deduced from the 
measured H c i(0). In principle, we evaluate the pen- 
etration depth A a b(0) using the expression: H t 



cl 



((/>o/47rA^ fc ) Ln/e, where k = A/£ is the Ginzburg-Landau 
parameter (here we assume that n is T-independent). Us- 
ing h£{°(0) = 695 Oe and k = 80, we obtain A q6 (0) ~ 
105 nm for sample No. 1 and 115 nm for sample No. 2. 
The value of A~ b 2 (0) , or equivalently the condensed carrier 
density n s /m* (superconducting carrier density/effective 
mass), allows us to check whether the well-known scal- 
ing behaviors between n s /m* and T c still works for the 
present system. In Fig. 4 we present our results to- 
gether with many others, including the FeAs-1111 sys- 
tem, cuprates, MgB2 and NbSe2 [3_7j. It is remarkable 
that Bao.6Ko.4Fe2As2 resides far away from the Uemura 
plot, which is contrasted by the case of the FeAs-1111 
system which is quite close to the Uemura plot. This dis- 
crepancy between FeAs-122 and FeAs-1111 implies very 
essential difference between the physics in the two sys- 
tems, this warrants further investigations. 

To summarize, we conduct magnetization measure- 
ments on Bao.6Ko.4Fe2As2 single crystals, and the lower 
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